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Abstract 

Complete next-to leading order QCD predictions for (2+1) jet cross 
sections and jet rates in deep inelastic scattering (DIS) based on a 
new parton level Monte Carlo program are presented. All relevant 
helicity contributions to the total cross section are included. Results 
on total jet cross sections as well as differential distributions in the 
basic kinematical variables x, W"^ and are shown for HERA en- 
ergies and for the fixed target experiment E665 at FERMILAB. We 
study the dependence on the choices of the renormalization scale hr 
and the factorization scale iip and show that the NLO results are 
much less sensitive to the variation of ^ = = than the LO 
results. The effect of an additional px cut to our jet definition scheme 
is investigated. 



I. INTRODUCTION 



The start-up of the HERA electron-proton colhder in 1992 marked the begin- 
ning of a new aera of experiments exploring Deep Inelastic Scattering of electrons 
and protons One of the topics to be studied at HERA will be the deep 

inelastic (= high Q^) production of multi jet eventsQ, where good event statistics 
are expected allowing for precision tests of QCD. Multijet events in DIS are first 
observed at the FERMILAB E665 experiment [Q. 

In this paper, we present complete results for (2+1) jet cross sections ("+" 
denotes the remnant jet) in DIS based on a new parton level Monte Carlo pro- 
gram DIS JET 0]. All helicity contributions to the total cross section (i.e. 
<7(7+l[(2 + 1) — jet] and crL[(2 + 1) — jet]; see below) are included. Note, that in 
contrast to jet production in e~^e~ experiments, it is not sufficient to calculate the 
contraction with the metrical tensor —g^^ on the hadronic tensor. For our kine- 
matical ranges at HERA (defined below), the ratio of ^^—g^i,H'^'^ / complete" is 

0. 75 (0.85) for the "low (high) range". Leading order (LO) and next-to leading 
order (NLO) matrix elements for the processes eP e+n jets+remant jet (n = 

1, 2, 3) are implemented in DIS JET. Using an invariant jet definition scheme in- 
troduced in we present results for total NLO (2+1) jet cross sections as well 
as differential distributions in the basic kinematical variables x, W"^ and Q^. The 
dependence on the renormalization scale hr and factorization scale yUi? is investi- 
gated. It is shown that the NLO predictions are much less sensitive to the choice 
of the scales that the LO results alone. Varying /i^ = = in the range of 
l/4p|, to 4 induces an uncertainty of roughly ±3% in the NLO (2+1) jet cross 
section predictions compared to about ±17% for the LO results for the "high 
Q'^ range" {Q"^ > lOOGeV^) at HERA. We have also studied the dependence on 
other choices of the scales [/i^ = O.IQ"^ — lOQ^; ycutW^] and found that these vari- 
ations lead to larger uncertainties (see table 1). However, a choice of the scales 
like /i^ ~ seems to be the more appropriate choice in the case of (2+1) jet 
production[|. We also explore the dependence of the (2+1) jet cross section on 
an additional cut on the transverse momentum px of the jets [pt is defined with 
respect to the 7* direction). It is shown that most of jets at a center of mass (cm) 



*Much higher effective luminosity will be achieved for low quasi-real photoproduc- 
tion of jets. The theoretical interpretation of these events, however, is more difficult, 
since there is added complication of the direct and the resolved photon contribution to 
the cross section, which are hard to separate 

^An extreme example would be (2+1) jet production at very low [Q"^ ~ 
IGeV^), where cannot be the relevant scale for high pT jet production. 
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energy of ^/s = 30 GeV (E665 experiment) are produced with a transverse mo- 
mentum less than 3 GeV using our jet definition scheme. Therefore, one should 
use resummation techniques to obtain a reliable perturbation expansion in this 
region. This problem is not tackled in this paper. Finally we also present results 
for different values of the jet resolution parameter i/cut- A comparison of jet mea- 
surements with our predicted QCD results provides a direct tool of determining 
as or Ams. Hadronization corrections may be minimized by restricting such an 
analyses to large momentum transfer which causes sufficiently large transverse 
momenta of the participating partons. 

The paper is organized as follows: In section II we discuss the general structure 
of the cross sections used in the calculation. Section III explains the jet definition 
scheme and in section IV we discuss kinematical ranges and numerical results for 
(2+1) jet cross sections and rates. Finally, section V contains a short summary. 

II. MATRIX ELEMENTS 

Consider deep inelastic electron proton scattering 

e^(/) + proton(P) proton remnant(j9r) + parton l{pi) . . . + parton n{pn) (1) 

Reaction proceeds via the exchange of an intermediate vector boson V = 
7*, Z, W. In this paper only the exchange by a virtual photon is considered. We 
denote the 7*-momentum by q, the absolute square by Q^, the center of mass 
energy by s, the square of the final hadronic mass by W'^ and introduce the 
scaling variables x and y: 

Q = 

s = 
W^ = 

X = 

y = 

At fixed s, only two variables in 

xW'^ = (1 - x)Q\ Q2 _ ^yg_ 

When one of the hadronic final state momenta pi is measured, reaction (|I|) is 
described in the one photon exchange by five parity conserved hadronic structure 
functions Hi — p 



—q^ = xys > 
(P + lf 

Pf = (P + g)2 (2) 



2Pq 



(0 < a; < 1) 



§ (0 < y < 1) 

(0) are independent, since e. g. 
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+ if4-^ {P'p'i + P^'P^) + i^s;^ (P^l^r - P'P^) (3) 

where we have introduced current conserved momenta variables = — ■^^^Q'^- 
The four /i <-> z/ symmetric structure functions i/i — contribute to so called 
T-even observables whereas the fi ^ u antisymmetric structure function H^, gives 
a contribution to "T-odd" observables p. To 0(as) in QCD one populates only 
Hi — H4 as there are no loop contributions to that order. LO contributions to 
these structure functions have been extensively studied in the literature ||5|,[7Hl4|. 



Note that in the totally inclusive case, where no final hadron momenta are 
measured, one only has contributions to Hi and H2 which are then denoted by 
the more familiar names Wi and W2- 

Experimentally one can measure the so called helicity cross sections <Ju+l, 
(Tt, 0"/, a A through lepton hadron correlation effects. They factorize the following 
y and (p dependence 

dcr[n - jet] ~ (1 + (1 - yf) dau+L[n - jet] - y^ daiin - jet] 

+ 2(1 — y) cos20 (icrT'[n — jet] — {2 — y)^\~ ycoscf) dai[n — jet] 

(4) 



+ yy^l -?/ sin rfcTA [n -jet] 

In eq. @ (p denotes the azimuthal angle between the parton plane (p, pi) and 
the lepton plane (/,/') (in the (7*-initial parton)-cms). 

The helicity cross sections ax {X E {U + L, L, T, I, A}) are linearly related 
to polarization density matrix elements of the virtual 7*. One has: 

<Ju+L ~ hoo + h++ + h— (5) 

ffL ~ /ioo (6) 

(Tt ~ + (7) 

o-/ ~ /i+o + /io+ " /^-o - /^o- (8) 

(^A ~ /i+o + ^0+ + ^-0 + ^0- (9) 

where hmm' = e*^{m)H'^^e^{m'), {m,m' = +,0,—) and e^(±)(e^(0)) are the 
transversal (longitudinal) polarization vectors of the 7* in the (7*— initial parton)- 
cms[]. Therefore (Tu+l and ai labels the unpolarized and longitudinal polar- 



■fXhe lepton hadron scattering process may be regarded as the scattering of a polarized 
off-shell gauge boson on the proton where the polarization of the gauge boson is tuned 
by the scattered lepton's momentum direction. 
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ization, cxr and cx/, a a mean transverse and transverse- longitudinal interference, 
respectively. These helicity cross sections are also linearly related to the five 
covariant structure functions defined in eq. (0) (see Apendix B of |^). 

In this paper we present results for (2+1) jet cross sections where we have 
integrated over the azimuthal angle (f). Therefore only (Ju^l and contribute 
in eq. (§). These two cross sections can technically be obtained by the following 
covariant projections on the (partonic) hadrontensor H^^ , which is calculated in 
fixed order perturbation theory (p = rjP denotes the momentum of the incoming 
parton and Xp = Q'^/{2pq)): 

/ 1 3x \ - 

au+L = l-^9f.u + -^P^^pA H^'in - jet] (10) 

= ^PmP^^ i^^ln - jet] (11) 
PI 

The 0(ag) (2+1) jet matrix elements represent a full NLO calculation includ- 
ing virtual and real corrections. The following subprocessess contribute to (2+1) 
jet cross sections up to 0(af): 

^^'^[tree, 0(a,) ] : 7* + g ^ g + G 
7* + G ^ g + g 
^^"[tree, 0(^2) ]:7* + g^g + G + G 

7* + g^g + g + g (12) 

7* + G^g + g + G' 
ii'''^ [virtual, 0(^2) ] : ^* + g g + G 

7* + G ^ g + g 

and the corresponding anti-quark processes with g g. Matrix elements for the 
complete contributions to (2+1) jet NLO O(a^) corrections are first discussed in 
p^ . The NLO matrix elements used in the MC program DISJET are based 
on these matrix elements. In addition, we have also included the full NLO scale 
dependent contributions. Note, that the second projection (~ V^iPv) m cr^j+i 
in eq. (p!0D gives a contribution of the order of 15 — 30% to the (2+1) jet cross 
sections depending on the kinematical ranges, wheras the contribution from ctl in 
eq. (|TID is fairly small (less than 1% in our kinematical ranges). This originates 
from the y dependent coefficients (1 + (1 — y)^) and —y^ and the fact that y 
is peaked at small values (see fig. 9 in 0]). The (2+1) jet NLO contributions 
originating from the projection with —g^l,v on the hadron tensor (see eq. ([T0|)) are 
first presented and discussed in detail in |]T6[. A complete list of tree level matrix 
elements with up to four partons in the final state can also be found in |T7| . 

The general structure of the NLO jet cross sections in DIS within the frame- 
work of perturbative QCD in given by: 
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da^-%2 + l)-]et] = JdT] faiv,fil) da^ip = r]P,aj{^^l),^il,^il) (13) 

where one sums over a = q,q,g. faijI^lA^) is the probabihty density to find a 
parton a with fraction r] in the proton if the proton is probed at a scale 
a"" denotes the partonic cross section from which colhnear initial state singular- 
ities have been factorized out at a scale /xl^ and implicitly included in the scale 
dependent parton densities fa{r], fi'^p). 

Let us briefiy discuss some technical matters that go into the NLO calcula- 
tion. The 0(ag) tree graph contributions in eq. ([T^p are integrated over the 
unresolved phase space region which are (2-1-1) jet like as defined in eq. (0). 
Infrared (IR) as well as colhnear (M) divergencies associated with the final state 
partons are cancelled against corresponding IR/M divergencies of the one-loop 
contributions. The remaining colhnear initial state divergencies are factorized 
into the bare parton densities introducing a factorization scale dependence ^p- 
Finally the ultraviolet (UV) divergencies are removed by MS renormalization 
which introduces a renormalization scale dependence /i/j. 



III. JET DEFINITION 

In order to calculate the (2-1-1) jet cross section we have to define what we call 
(2-1-1) jets by introducing a resolution criterion. As has been elaborated in detail 
in 1^,1^ energy-angle cuts are not suitable for an asymmetric machine with its 
strong boosts from the hadronic cms to the laboratory frame. As a jet resolution 
criterion we use the invariant mass cut criterion defined in [^^|15|,|16| such that 



Sij = {Pi + Pjf > = max [vcutMl, M^^ (z, j = 1, . . . , n, r; z 7^ j) 

(14) 

where Ucut is the resolution parameter and Sij is the invariant mass of any two 
final state partons, including the remnant jet with momentum pr = {1 — v)P- 
Mc is a typical mass scale of the process which defines the jet definition scheme. 
In this paper we choose = W^. This corresponds to the "PF-" scheme in [^. 



Other jet definition schemes based on a kx algorithm are proposed in |]T8|. Mq 
is a fixed mass cut which we have introduced in order to clearly separate the 
perturbative and non-perturbative regime in the case where is small. Mq is 
fixed to 2 GeV H7| in all our results. 
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IV. NUMERICAL RESULTS 



We will now turn to our numerical cross section results and present results 
for the actual HERA cm energy of 295 GeV as well as for the FERMILAB fixed 
target experiment E665 with a cm energy of 30 GeV. All results are based on a 
new Monte Carlo program DISJET The Monte-Carlo routines are using the 
VEGAS-package for numerical integration. Parton distributions are incorpo- 
rated from the packages p(]| , pT| . Our standard set of parton distribution functions 
is MRS set D- [^. If not stated otherwise, we use the one-loop formula for the 
strong coupling constant for our LO results 



R) 



(33 - 2nf) In ^ 
whereas we employ the two loop formula 



(15) 



a 



sMS\t^R 



127r 



(33-2n;)ln(/iyA2 



1 - 



6(153- 19n/) In ln(/iyA2 
(33-2n^)2 ln(/iyA2) 



(16) 



in our NLO predictions. The A^g value is chosen consistent to the value from 
the parton distribution functions. The value of is matched at the thresholds 
q = rriq and the number of flavours in is fixed by the number of flavours 
for which the masses are less than /i^. Furthermore the number of quark flavours 
that can be pair-produced are set equal to n/ chosen in as- In all predictions, the 
renormalization scale and the factorization scale are set to be equal: /i^ = = 
/i^. For the flne structure constant a we adopt the running coupling formula. 
The following kinematical cuts are used for the HERA results: 



0.001 
0.04 

600 GeV^ 



< 
< 
< 



X 

y 



< 
< 



0.1 
0.95 



(17) 



In addition, we use two different ranges: 
"Low g2 range:" 4 GeV^ < < 100 GeV^ 
"high Q2 range:" 100 GeV^ < 

Note that the kinematical cuts in eq. (0) are not independent, for example the 
X and y cut also imply W > 55 GeV (W"^ = (1 — x)ys) and Xmax in the low 
range is Xmax = 0.0287 rather than 0.1 (Q^ = xys). 

Our kinematical range for the E665 experiment is deflned by 0]: 



0.003 
0.08 
4 GeV^ 
400 GeV^ 



< 
< 
< 
< 



X 

y 



< 
< 



0.95 
25 GeV^ 
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In fig. 1 we show the dependence of the (2+1) jet rate on the resolution 
parameter Ucut for the two different ranges at HERA (a,b) and for the E665 
experiment (c). Solid (dashed) lines correspond to NLO (LO) predictions. The 
renormalization scale and the factorization scales are set equal to (lower 
curves) and (upper curves). The transverse momentum pt of the jets is defined 
in eq. (|l^). One observes, that the NLO corrections lower the (2+1) jet rate by 
about 15%. At Vcut = 0.02 the (2+1) jet rate exceeds about 4% (12%) for the low 
(high) range at HERA. For the E665 experiment, the (2+1) jet rate exceeds 
about 10-11% for i/cut = 0.04. In the following, we use ycut = 0.02 for HERA and 
Vcut = 0.04 for E665 as our standard values. 

To get a feeling for the theoretical uncertainties originating from the choice 
of the renormalization and factorization scales, we show numerical results for the 
(2+1) jet cross sections for different fj,"^ values in table 1. Let us first comment on 
the choice of the scales. In DIS scattering it is natural to take /x^ = /x^ = /i^^ = Q^. 
However, for jet production, the transverse momentum p^ of the jets should also 
be considered as a relevant scale. The transverse momentum is defined with 
respect to the 7* direction. For the (1+1) case, one has two jets, the remant jet 
and the struck parton jet, both with zero pj^. In our case of (2+1) jet production, 
one expects two partonic jets with nearly opposite large px and the remnant jet 
at Pt = 0. For the LO processes, px is given by: 

p't = Q'^—^z{l - z) (19) 

Xp 

with 

xp = ^ = ^- . = ^ (20) 
2pq 7] pq 

where pi is the four momentum of one of the final partons. Finally, in analogy 
to the jet analysis in e^e~ experiments one may use HcutW^ (see. eq. (p!^)) as 
a possible scale. To avoid to small scales for perturbation theory, the scales are 
clipped at a minimum value of 2 GeV^. Varying = = between Q^, 
0.25 — Ap"^ and UcutW^ induces an uncertainty of roughly ± 17% (± 4%) in 
LO (NLO) for the high range at HERA and ± 20% (± 10%) in LO (NLO) for 
the E665 experiment. Therefore, the uncertainty in the theoretical predictions is 
markedly reduced by the NLO corrections. 

In figs. 2-10 we show the dependences of the total cross section and the (2+1) 
jet cross sections as well as the (2+1) jet rate on the basic kinematical variables 
X, and . Figs, a (b) are for the low (high) Q"^ range at HERA and c shows 



The intrinsic transverse momentum of the partons in the target is neglected 
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predictions for the E665 experiment. LO results are given in figs. 2,3,6,7,11-15 
whereas figs. 4,5,8,9 and 10 show NLO predictions. 

In order to estimate the theoretical uncertainty from the choice of the scales, 
all results are given for jj? = n\ = ij?p = ap"^ (a=l/4, 1,4: solid lines), /i^ = 
aQ^ (a=l/10, 1, 10: dashed lines) and fj? = VcutW^ (dotted lines). 

Let us first comment on the W distributions shown in figs. 2-5. The (2+1) jet 
rate decreases with increasing W . This is mainly a reflection of our jet definition 
in eq. (|I^) where the required invariant mass of the jets increases with increasing 
W'^. The (2+1) jet rate for the high Q"^ sample is larger than for the low Q'^ 
range since the larger causes at the average larger transverse momenta of the 
participating partons (see eq. (p!9D) and larger invariant masses. Comparing figs. 
2,3 and figs. 4,5 one observes again that the uncertainty from the choice of the 
scales in the LO predictions is very large (figs. 2,3) whereas this uncertainty is 
markedly reduced by including the NLO corrections (figs 4,5). 

Turning now to the x distributions one observes a quite different behaviour 
of the cross sections for the different kinematical ranges in figs. 6 and 8. The 
X distributions are mainly governed by the behaviour of the respective parton 
distributions in the allowed kinematical regions. Note that Xmax = 0.0287 for the 
low range (figs 6a and 8a). The (2+1) jet rates in figs. 7 and 9 are increasing 
with increasing x. This is mainly an effect of the increasing (2+1) jet phase 
space in our jet definition scheme (for a detailed discussion of the 0(0;^) (2+1) 
jet phase space see 0). Comparing the distributions in figs. 6-9, one observes 
again that the NLO predictions in figs. 8,9 are much less scale dependent than 
the LO predictions in figs. 6,7. 

Fig. 10 shows NLO predictions for the (2+1) jet rate as a function of ^/Q^. 
Results are given for QsifJ'^ = Q^) (lower solid curve) asifi'ji = Pt), (upper solid 
curve) and «<, = fix = 0.25 (dotted curve). Note that the (2+1) jet rate is an 
increasing function with increasing Q^. Therefore the or dependence in 
as is overcompensated by the increasing (2+1) jet phase space relative to the 
total cross section in our jet definition scheme. One observes sizable differences 
between the predictions using a constant and scale dependent (two- loop) a^. 
Therefore a clear discrimination between the solid and dotted curves should be 
possible with the expected event statistics at HERA. 

Note also, that /i^ = p"^ tends to predict larger (2+1) jet rates and differential 
distributions than //^ = (see also table 1 and figs. 1-10). This is in particular 
true for lower W'^ and higher x values. 

In order to explore the pt dependence of the jet rates in more detail, fig. 11 
shows results for the (2+1) jet rate as function of an additional pt cut. Applying 
a Pt cut of 4 GeV reduces the (2+1) jet rate for the low (high) range at 
HERA from ~ 4% to 2-3% (~ 12% to ~ 10%) whereas the (2+1) jet rate falls 
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from 10% to less than 1% for the FERMILAB experiment. Fig. 11c shows 
also, that more than 50% of the jets at ^/s = 30 GeV are produced with a 
transversal momentum less than 2 GeV. However, these values are too small to 
allow for reliable predictions in fixed order perturbation theory and resummation 
techniques should be used in this kinematical region. In fact , the experimental 
results presented by the E665 collaboration are significantly higher than our 
second order predictions in figs. 1-9 c. 

In table 2, we give results for jet cross sections with an additional prmin cut. 
Figs. 12-15 show the corresponding W and x distributions. Note that the (2 + 1) 
jet rates in figs. 12-15 are less x and W dependent than our results without the 
additional PTmm cut. This demonstrates that the increase of the jet rates for 
small W (figs. 3,5) and large x values (figs. 7,9) is mainly an effect of low px 
jets. This is also clear by an inspection of eq. ([T9|) . 

V. SUMMARY 

The (2+1) jet production in DIS is calculated up to NLO in perturbative QCD 
using an invariant jet definition scheme. The theoretical uncertainties originating 
from variations of the renormalization/factorization scales are well under control. 
Differential distributions of (2+1) jet rates in W'^,x and are presented. A 
comparison with the large numer of events expected at HERA in the near future 
will allow for precision tests of perturbative QCD. It is shown that an additional 
Pt cut to our jet definition scheme is necessary to obtain reliable predictions for 
jet production at the energy of the E665 experiment at FERMILAB. 
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TABLES 





low Q^: LO 


low Q^: NLO 


high Q2. lo 


high NLO 


E665: LO 


E665: NLO 




2360 


1930 


600 


547 


1200 


1020 


10 • 


1980 


1930 


472 


502 


879 


932 


0.1 • Q2 


3050 


2000 


770 


548 


1660 


1130 




2080 


1960 


602 


553 


1040 


983 




2320 


2020 


719 


579 


1310 


1120 




2100 


2000 


615 


562 


1040 


1010 


0.25 • 


2670 


2140 


857 


583 


1590 


1180 



TABLE L (2+1) jet cross sections in [pb] for different choices of the renormalization 
and factorization scales (column 1). LO and NLO results are shown for the "low Q^" 
(column 2 and 3) and the "high Q^" (column 4 and 5) ranges at HERA {ycut = 0.02) 
and for the E665 experiment (column 6 and 7) (ycut = 0.04). 





low Q2 


high Q2 


E665 


Pt 


1140 


451 


75 


^■pI 


1040 


394 


59 


0.25 • 


1230 


514 


103 



TABLE IL LO predictions for (2+1) jet cross sections in [pb] with an additional 

PTmin cut of 4 GeV. The two-loop formula for Ug is used in the calculation. 
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FIGURE CAPTIONS 



Fig. 1. (2+1) jet fraction versus the cut variable Ucut for the low and high Q"^ range 
at HERA (a,b) and for the E665 experiment (c). The kinematical ranges 
are defined in eqs. (0,|1^). Solid (dashed) lines correspond to NLO (LO) 
predictions. In the upper (lower) lines the renormalization scale and 
factorization scale hf are set equal to /i^ = {^"^ = Q"^). 

Fig. 2. W dependence of the total cross section (upper line) and the (2+1) jet cross 
section in LO for the low (a) and high (b) range at HERA and for 
the E665 experiment (c). ycut = 0.02 (0.04) in a,b (c). The different curves 
for the jet cross section belong to different choices of the renormalization 
scale fiR and factorization scale fip: 

solid lines: /i^ = = = l/4p|,, p"^, Ap\ (from top to bottom), 
dashed lines: /i^ = 0.1(5^,<5^, lOQ^ (from top to bottom), 
dotted line: /i^ = VcutW^. 

Fig. 3. W dependence of the (2+1) jet rate in LO. Solid, dashed and dotted lines 
as in fig. 2. 

Fig. 4. same as fig. 2, but NLO predictions. The two upper lines in c) show the 
LO (dashed) and NLO (0(q;s), solid) predictions for the total cross section. 
The difference between these results for HERA energies are too small to be 
visible in a) and b). 

Fig. 5. same as fig. 3, but NLO predictions. 

Fig. 6. X dependence of the total cross section (upper line) and the (2+1) jet cross 
section in LO for the low (a) and high (b) range at HERA and for 
the E665 experiment (c). Solid, dashed and dotted curves as in fig. 2. 

Fig. 7. X dependence of the (2+1) jet rate in LO. Solid, dashed and dotted curves 
as in fig. 6. 

Fig. 8. same as fig. 6, but NLO predictions. 

Fig. 9. same as fig. 7, but NLO predictions. 

Fig. 10. y/Q^ dependence of the (2+1) jet rate in NLO for [j? = (lower solid line) 
/i^ = p'^ (upper solid line). Also shown is the result for = const. = 0.25 
(dotted line). 

Fig. 11. LO (2+1) jet rate as a function of prmin for the low and high Q"^ range at 
HERA with i/cut = 0.02 (a,b) and for E665 experiment (c) {i/cut = 0.04). 
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The three curves are for fj,'^ — /jj^ = fxj^ = Ap^ (from top to 

bottom). We have used the two- loop formula for Ug. 

Fig. 12 W dependence of the total cross section (upper line) and the (2+1) jet cross 
section in LO for the low (a) and high (b) range at HERA and for 
the E665 experiment (c) with an additional PTmin cut of 4 GeV (for the 
jet cross sections). Same parameters as in fig. 11 for the (2+1) jet cross 
sections. 

Fig. 13 W dependence of the (2+1) jet rate in LO with an PTmin cut of 4 GeV 
(parameters as in fig. 11). 

Fig. 14. X dependence of the total cross section (upper line) and the (2+1) jet cross 
section in LO for the low (a) and high (b) range at HERA and for 
the E665 experiment(c) (parameters as in fig. 11). 

Fig. 15. X dependence of the (2+1) jet rate in LO with an PTmin cut of 4 GeV 
(Parameters as in fig. 11). 
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